Superradiant Phase Transitions and the Standard Description of Circuit QED 
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We investigate the equilibrium behaviour of a superconducting circuit QED system containing 
a large number of artificial atoms. It is shown that the currently accepted standard description 
of circuit QED via an effective model fails in an important aspect: it predicts the possibility of 
a superradiant quantum phase transition, even though a full microscopic treatment reveals that a 
no-go theorem for such phase transitions known from cavity QED applies to circuit QED systems 
as well. We generalize the no-go theorem to the case of (artificial) atoms with many energy levels 
and thus make it more applicable for realistic cavity or circuit QED systems. 
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Recent years have seen rapid progress in fabrication 
and experimental control of superconducting circuit QED 
systems, in which a steadily increasing number of artifi- 
cial atoms interact with microwaves [EE]- These devel- 
opments set the stage to study collective phenomena in 
circuit QED. An interesting question in that context is 
whether a system with many artificial atoms undergoes 
an equilibrium quantum phase transition as the coupling 
of artificial atoms and electromagnetic field is increased. 
Quantum phase transitions of this type have been in- 
tensely discussed for cavity QED systems [5-10J and 
are known as superradiant phase transitions (SPTs) [6]. 
However, in cavity QED systems with electric dipole cou- 
pling their existence is doubted due to a no-go theo- 
rem [8 J . Recently, it has been claimed that SPTs are pos- 
sible in the closely related circuit QED systems with ca- 
pacitive coupling jlQHT2] . This would imply that the no- 
go theorem of cavity QED does not apply and challenges 
the well-established analogy of circuit and cavity QED. 

Here, we show in a full microscopic analysis that circuit 
QED systems are also subject to the no-go theorem. We 
argue that such an analysis is necessary since the stan- 
dard description of circuit QED systems by an effective 
model (EM) is deficient in the regime considered here. 
A toy model is used to illustrate this failure of an EM. 
Finally, we close a possible loophole of the no-go theorem 
by generalizing it from two-level to multi-level (artificial) 
atoms. Thus, our work restores the analogy of circuit and 
cavity QED and rules out SPTs in these systems under 
realistic conditions that have not been covered before. 

The Dicke Hamiltonian in cavity and circuit QED. — 
Both circuit QED systems and cavity QED systems with 
N (artificial) atoms (Fig. [I]) are often described by the 
Dicke Hamiltonian [13J (h = 1) 
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with energy splitting Q between ground state \g)k = ( ? ) k 





FIG. 1: Cavity QED system with N atoms (a) and circuit 
QED system with N Cooper-pair boxes as artificial atoms (b) . 

and excited state \e)k = (^^2 are Pauli matri- 

ces). In the case of circuit QED, we assume Cooper-pair 
boxes as artificial atoms, which justifies the two- level ap- 
proxiation. Our main results, though, hold for arbitrary 
charge-based artificial atoms (capacitive coupling) [T4] . 
Further, a) generates a photon of energy u. Matter and 
field couple with a strength A. The ft-term, often ne- 
glected in other contexts, will become crucial below. In 
cavity QED, Hr> is obtained from minimal coupling of 
atoms and electromagnetic field. For a single atom (with 
n electrons) at a fixed position, the Hamiltonian reads 
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+ V int (ri,...,r n ). (1) 



The pA- and A 2 -terms in the analog TV-atom Hamilto- 
nian yield the A- and ft-term in Hd, respectively. In cir- 
cuit QED, Hb arises from a widely used EM for a charge- 
based artificial atom in a transmission line resonator [15] , 

n° ciI = 4£c5>-*>)»H - § £(|^+1>H+H.c). 

Here, v counts the excess Cooper pairs on the island, 
Ej and Eq = e 2 /(2(Cq + Cj)) are the Josephson energy 
and the charging energy of the Cooper-pair box, Cq and 
Cj are the coupling capacitance and the capacitance of 
the Josephson junction. Moreover, v — Cq(Vq + V)/2e, 
where Vq is an external gate voltage and V the quan- 
tum voltage due to the electromagnetic field in the res- 
onator. The Cooper-pair box is assumed to be operated 
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at its degeneracy point [15]. As it is described by macro- 
scopic quantities (like Eq) and only a single degree of 
freedom (V), 1-Uq t is an EM for a Cooper-pair box in a 
transmission line. Starting either from H|? av or from H[? ir , 
one obtains Hb using the following approximations: The 
N (artificial) atoms are treated as identical, noninter- 
acting two- level systems with ground states \g) and ex- 
cited states |e) which are strongly localized compared to 
the wavelength of the single considered field mode (i.e., 
A(if ) * A = A)e(at + a) and V(r k ) « V = V (at + a)). 

Superradiant phase transitions and no-go theorem. - 
In the limit TV — >> oo, undergoes a second order quan- 
tum phase transition at a critical coupling strength [SHE] 
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This phase transition was discovered by Hepp and Lieb 
for with k = and termed SPT [6]; see [9] for re- 
cent studies. At A c , the atoms polarize spontaneously, 
(Efc o-JJ/iV ^ -1, and a finite macroscopic photon occu- 
pation arises , (a^a)/N ^ 0. This quantum critical point 
is signaled by a gapless excitation (Fig. [2ja)). 

In cavity QED systems, however, A c cannot be reached 
if the ft-term is not neglected |8J . That is because A and 
k are not independent of each other. Let us define a 
parameter a via n = aX 2 /Q. Then Eq. Q becomes 
A 2 (l — a) = and criticality requires a < 1. With 

Aq — l/y / 2eocJF (V is the volume of the cavity) one finds 
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where d 

But the Thomas-Reiche-Kuhn sum rule (TRK) 
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for the Hamiltonian H° = EILiPf /2wi+Vint(ri, . . . , r n ) 
of an uncoupled atom with spectrum |Z)} implies 
Q\e • d| 2 < ne 2 /2m, consequently a cav > 1. This is 
known as the no-go theorem for SPTs [HI US]. Notice 
that a cav determines how strongly Q\e • d| 2 exhausts the 
TRK. We remark that a direct dipole-dipole coupling be- 
tween atoms (omitted here) can lead to a ferroelectric 
phase transition, which however occurs only at very high 
atomic densities |16| . 

Surprisingly, the no-go theorem was recently argued 
not to apply in circuit QED [10]. Indeed, the standard 
EM of circuit QED yields 
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where L denotes the length of the transmission line res- 
onator, c its capacitance per unit length, and we have 
used V = {yjjLcfl 2 [15J. Here a cir = Ej/4E C < 1 is 
easily possible [lj. According to this argument, a SPT 
should be observable in a circuit QED system. 




FIG. 2: (a) Excitation energies e+ and e_ of the Dicke 
Hamiltonian Tin versus coupling A (in units of a; = Q), for 
a = kQ/X 2 = 0, 0.8, 1, 1.2. For a = 0, e_ vanishes at A = 0.5, 
thus signaling a SPT. Only a > 1 is compatible with the TRK 
sum rule. For these a, e_ — > y/l — 1/a and remains finite for 
all A. The excitations e±(A) of 7it m correspond to a = 1. 
(b) Toy model of an (artificial) atom. The blue line indicates 
the degree of freedom in the simplified effective model. 

Effective models and superradiant phase transitions. — 
The EM has proved to be a very successful description 
of circuit QED whose predictions have been confirmed in 
numerous experiments. However, the circuit QED setups 
operated so far contained only few artificial atoms. It is 
not obvious that an EM provides also a good description 
of circuit QED systems with N ^> 1 atoms and, thus, a 
proper starting point to study SPTs in circuit QED. We 
now present a toy model illustrating how an EM similar 
to the one in circuit QED erroneously predicts a SPT. 

The toy model consists of N harmonic oscillator po- 
tentials with frequency each trapping n noninteract- 
ing fermions of mass m and charge e, which all couple to 
a bosonic mode with frequency uu (Fig. [2jb)). This toy 
model can be viewed as a very simplified description of 
(artificial) atoms with n microscopic constituents inside 
a resonator. It is governed by the Hamiltonian 
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where we assume again A(x k ) w A = Aq{o)-\-o). Since A 
couples only to the center of mass coordinate of the kth. 
harmonic oscillator, T~L can be diagonalized exactly (|17J): 

nN—l 

n tm = e ± (ala ± + -) + ^ + -), (7) 



2e^(A) = cj 2 + 4.KUJ + n 2 ± vG^T ^uj-n 2 ) 2 + 16X 2 ujn. 

Here, generate excitations that mix photon field with 
collective center of mass motion, the b\ excite the remain- 
ing degrees of freedom, A = AoVtd\fN and n = X 2 /ft. As 
d = (n\ex\n — 1) = e^/n/2mQ^ the TRK is exhausted. 
Note that e±(A) are also the relevant excitation energies 
of for N — > oo, as can be shown using methods of 
Ref. [9] (see [17]), and demanding e_ = yields Eq. 
One sees that e±(A) is real and nonzero for all A and that 
the ground state energy is an analytic function of A (Fig. 
|2^a)). Hence, no quantum phase transition is possible. 
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Let us now consider an EM for the toy model. Sim- 
ilar to the standard EM of circuit QED, we focus on 
the fermion with the highest energy in the kth har- 
monic oscillator and treat it as two-level system with 
\g k ) = \n- l) k and \e k ) = \n) k (Fig. [2^b)). Accounting 
only for one fermion per 'atom' yields a Dicke Hamil- 
tonian with A em = A and ft em = n/n = Ne 2 Aq/2tu. 
Crucially, only A em depends on n. This allows A em to be 
increased at constant ft em , therefore a em = 1/n can be 
< 1 and a SPT is possible. This failure of the EM can be 
interpreted as follows. The relation A = A em oc d oc y/n 
reveals that the coupling of an 'atom' to the bosonic mode 
is fully captured by the EM and grows with atom size n. 
However, in a proper description of the system, increasing 
the coupling by increasing n unavoidably also increases 
n in proportion to n: all fermions of all atoms couple to 
the bosonic mode and each causes an A 2 -term. This is 
lost in the EM with only one degree of freedom per atom. 
Interestingly, a em < 1 only if n > 1, i.e., as long as the 
effective description actually neglects degrees of freedom. 

Microscopic description of circuit QED. — This exam- 
ple suggests not to rely on the standard description for 
investigating SPTs in circuit QED. Although the dipole 
coupling of field and qubit states might be fully repre- 
sented by A C i r , ft C i r could still underestimate the A 2 -terms 
of all charged particles in the Cooper-pair boxes. Instead, 
let us describe a circuit QED system with N artificial 
atoms by a minimal-coupling Hamiltonian that takes all 
microscopic degrees of freedom into account: 
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As we allow arbitrary charges q k and masses m\ and 
an arbitrary interaction potential Vj nt of the n k con- 
stituents of the kth artificial atom, H m ic captures the 
coupling of N arbitrary (but mutually noninteracting) 
objects to the electromagnetic field in the most general 
way. We subject it to the same approximations that led 
from H° ir , the EM of circuit QED, to For iden- 

tical artificial atoms {n k , q k ,m k } — > {n,^,mj. The 
Hamiltonian of an uncoupled artificial atom then reads 
#mic = Pf /2m, + y int (r l5 . . . , r n ). Notice that the 

qubit states \g) and |e) of the artificial atom, which in the 
standard EM are superpositions of the charge states |^), 
are among the eigenstates {|/)} of ^mic By expanding 
Hmic in the {\g)k, |e)fe}-basis and taking A(r^) « A, one 
obtains the Dicke Hamiltonian %d with 
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and d = (e\ J27=i ^ r d#)- This microscopic description of 
circuit QED facilitates the same line of argument which 
in Ref. |8 allowed the conclusion that there is no SPT 
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in cavity QED: Criticality (Eq. (|2|) requires Q\e • d| 
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Hence, the no-go theorem of cavity QED applies to circuit 
QED as well. This result confirms the analogy of cavity 
and circuit QED also with respect to SPTs. It has been 
obtained under the same approximations that led from 
the standard description of circuit QED, H^ r , to Hb with 
A c i r and ^ c i r . The discrepancy of the predictions of the 
microscopic and the standard description of circuit QED 
thus shows the limitations of the validity of the latter. 
This might be important for future circuit QED archi- 
tectures with many artificial atoms in general, even for 
applications not related to SPTs. We emphasize, though, 
that our conclusion neither forbids SPTs in circuit QED 
systems with inductively coupling flux qubits [19], nor is 
it at odds with the great success of the standard descrip- 
tion for few-atom systems: there, the deficiency of n c i v 
does not manifest itself qualitatively as the ft-term in Hb 
mimics slightly renormalized system parameters Co and A. 

Possible loophole in the no-go theorem. — Although the 
two-level approximation for the anharmonic spectrum of 
(artificial) atoms is well justified in many cases, one might 
argue that higher levels should be taken into account in 
this context. Indeed, a SPT does not require Q w co>, and 
thereby does not single out a particular atomic transition. 

For a more profound reason for dropping the two- 
level assumption, consider the elementary question of 
how a resonator's frequency uj is shifted by the pres- 
ence of N mutually noninteracting atoms. This sit- 
uation is described by H m ic It can be rewritten as 
U mic = wata + £f=i {HK C + U k pA + n%), where H k pA 
and l~L k A 2 are the pA- and the A 2 -terms due to the kth. 
atom ([U]). Let us perturbatively calculate the fre- 
quency shift 5uo = Su p A + Sujaz caused by both ^ 
and (HZ]). To this end, take uj <C Vt k m for all m, k, 

where is the rath excitation energy of H^ ic . Remark- 
ably, it turns out ([II]) that Suj p a (< 0) and 8uja 2 (> 0) 
cancel almost exactly by virtue of the TRK, and the to- 
tal frequency shift 5uo is small, 5uo ~ (uu/Q^) 2 . As a 
SPT corresponds to Suj = —a;, the significance of both 
pA- and A 2 -terms for the existence of a SPT becomes 
clear. The pA-terms cause a strong negative frequency 
shift and favor a SPT, the A 2 -terms do the opposite. 
This means, most crucially, that one must not truncate 
pA- and A 2 -terms in an unbalanced way for assessing 
the possibility of a SPT by an approximate Hamiltonian. 
Dropping the A 2 -terms in Hb (^ = 0) leads to the pre- 
diction of a SPT. In contrast, %b with nonzero n fully 
incorporates the A 2 -terms of H m ic- But, due to the two- 
level approximation, it allows for only one matrix element 
of the pA-terms per atom, thereby possibly underesti- 
mating the tendency towards a SPT. To exclude SPTs 
in cavity and circuit QED, a generalization of the no-go 
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FIG. 3: Situation of the generalized no-go theorem. Many 
multi-level (artificial) atoms couple to the photon field. Tran- 
sitions between excited atomic states are irrelevant for the 
low-energy spectrum of the system. 

theorem to (artificial) atoms with more than two energy 
levels is necessary. 

Generalized no-go theorem. — Let us consider N — > oo 
identical atoms couped to a field mode with frequency uo. 
The atomic Hamiltonians H^ ic ma y nave an arbitrary 
spectrum {fij, \l^) = \l)k}, with = and \i excited 
states (Fig. With d hV = e • (l\ £? =1 ^|Z'), the full 
Hamiltonian of the system reads 



H m ic vat a + n(o) + a) 2 + ^ ^ (^i^i, v\h){h\ 

k=l 1,1' =o 

+ iAj(fii/-ftiKi/(ot + o)|! fc )(jj;i). (io) 



We now follow a strategy similar to that of Refs. [9J: We 
derive a generalized Dicke Hamiltonian HgD having the 
same low-energy spectrum as H m i c for a small density of 
atoms, N/V ~ 0, using Aq oc V -1 / 2 as small parameter. 
We then check whether H g D has a gapless excitation if 
the density is increased, which would signal a SPT and 
mark the breakdown of the analogy of H g D and H m i c . 

Expanding the eigenstates and -energies of T^mic &s 
\£) « E^o^ol^) and £ oc £ ss , ^o +S '(^mic|^>, we 
note that contributions from all d^o,z'/o terms may be 
neglected: they are smaller than those retained by a fac- 
tor of at least one power of Aq (for s-\-s f > l)or£/7V<C 1 
(for s + s' < 1), where f = ^ k ^2 l>0 l(^l^o)| 2 is the num- 
ber of atomic excitations in |£q), which is <C N for low- 
lying eigenstates ( [1 TJ ) . We thus define H g D by setting 
d^o^'/o ^ in Hmi C . Up to a constant, we find ([T7]) 



^n / ^ + ^^(^ t +&/)(a t +a), (11) 
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by introducing b\ = ^= J2k=i l^)(0fc| as collective exci- 
tation, omitting the energy of the 'dark' collective exci- 
tations ([17]), and removing the ft-term by a Bogolyubov 
transformation yielding uj — » a; = Vo; 2 + and — >• 
A/ = i/fA/, with X t = A Qi\d ,i\y/N. In the limit of 
dilute excitations, the bi are bosonic, [6/ , 6 j", ] = |18| . 
The system undergoes a SPT if an eigenfrequency of 
H g D can be pushed to zero by increasing the couplings 
A;. We cannot calculate the e^'s explicitly, but we will 
show that the assumption e$ = contradicts the TRK. 
An €i solves the characteristic equation (|17|) 
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and contradict the TRK for H^ ic (Eq. ([9|), which ensures 
that the right side is negative even if the entire atomic 
spectrum is incorporated. This result is irrespective of 
the details of the atomic spectra. Note that for n = 0, 
the negative term on the right side of Eq. ( |A.46 ) vanishes, 
and one recovers the SPT for critical couplings A/ c with 
JXi A ?cM = ^/ 4 - This resembles Eq. ^ with k = 0. 

Experimental evidence for our conclusions could be 
gained by simply probing the shifted resonator frequency 
of a suitable circuit QED system. Consider a sample con- 
taining N artificial atoms with X/\fN = 2tt x 120MHz 
and n/2ir = uj/2tt = 3GHz. If a cir = E 3 /4E C = 0.1, as 
predicted by the standard theory of circuit QED, the sys- 
tem should exhibit signatures of criticality for N = 174 
(according to Eq. Q), and its resonator frequency should 
be close to zero. But even if we assume a = 1, the mini- 
mal value compatible with the TRK (that corresponds to 
ideal two- level atoms), we find the lowest excitation e_ to 
be still at e_ ~ 2ir x 2GHz. We expect these phenomena 
to be insensitive to small fluctuations of the atomic pa- 
rameters [19J and hence to be experimentally observable. 

We thank S. M. Girvin, A. Wallraff, J. Fink, A. Blais, 
J. Siewert, D. Esteve, J. Keeling, P. Nataf, and C. Ciuti 
for discussions. Support by NIM, the Emmy-Noether 
program, and the SFB 631 of the DFG is gratefully ac- 
knowledged. 



[I] A. Wallraff et al, Nature 431, 162 (2004). 

[2] J. M. Fink et al, Phys. Rev. Lett. 103, 083601 (2009). 

[3] M. Neeley et al, Nature 467, 570 (2010). 

[4] L. DiCarlo et al, Nature 467, 574 (2010). 

[5] J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 

Phys. 73, 565 (2001). 
[6] K. Hepp and E. H. Lieb, Ann. Phys. 76, 360 (1973). 
[7] Y. K. Wang and F. T. Hioe, Phys. Rev. A 7, 831 (1973). 
[8] K. Rzazewski and K. Wodkiewicz, Phys. Rev. Lett. 35, 

432 (1975). 

[9] C. Emary and T. Brandes, Phys. Rev. Lett. 90, 044101 
(2003); C. Emary and T. Brandes, Phys. Rev. E 67, 
066203 (2003). 

[10] P. Nataf and C. Ciuti, Nature Commun. 1, 72 (2010). 

[II] G. Chen, Z. Chen, and J. Liang, Phys. Rev. A 76, 055803 
(2007). 

[12] N. Lambert et a/., Phys. Rev. B 80, 165308 (2009). 
[13] R. H. Dicke, Phys. Rev. 93, 99 (1954). 
[14] J. Clarke and F. Wilhelm, Nature 453, 1031 (2008). 
[15] A. Blais et a/., Phys. Rev. A 62, 062320 (2004). 
[16] J. Keeling, J. Phys.: Condens. Matter 19, 295213 (2007). 
[17] For details, see Appendix. 
[18] J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 
[19] P. Nataf and C. Ciuti, Phys. Rev. Lett. 104, 023601 
(2010). 



5 



Appendix 



We provide intermediate steps for the derivation of some important statements and equations 
of the main text. For clarity, formulas contained in the main text are typeset in blue. 



Thomas-Reiche-Kuhn sum rule. We derive the TRK [AlJ for the Hamiltonian 

n 2 



i=l 



yielding Eq. (|9| of the main text; Eq. Q follows as a special case. The derivation of the TRK 
is based upon the identities 

n 2 n n 71 71 

£ = A e ■ £ q ^ ' e • £ f^l ' £ I? = * £ * r *l ' (A - 15) 

2=1 2=1 i' = l 2=1 2=1 

for a real unit vector e. We denote the eigenspectrum of H^ ic by {E^ \l)}. It comprises a 



ground state \g) of energy E g . The TRK follows by combining the commutators of Eqs. (A.15): 

n 2 n n 

£ ^7 = <9\ [e ■ E ' I • Kic E «H] l5> ( A ' 16a ) 

2 = 1 % 2 = 1 2 V = 1 

n 

= -^Msl^r^ 2 . (A.16b) 



Diagonalization of and 'Htm* ^ is demonstrated that the diagonalization of both 
the Dicke Hamiltonain Hb for N — >> oo and the Hamiltonian Htm describing the toy model 
can be reduced to the diagonalization of special cases of H g D, which appears in the context of 
the generalized no-go theorem. The characteristic equation of H g D, that will be derived in the 
last section of this appendix, is solvable for the special cases and yields the diagonal forms of 
Kb and H tm . 



7-Lb= ua)a + — ^ a 



Diagonalization of First, we focus on the Dicke Hamiltonian 

^ N a AT 



: uja^a + 



■£*S 

k=l 



^cr^(a f +a) + ^(a + a t ) 2 

/e=l 
iV 

5^(7*(a+ + a) + C 



fc=l 



(A.17a) 
(A.17b) 



with Co = \Juo 2 + 4^cj, A = y/oo/tiX, and C = (a) — gj)/2. The Hamiltonian (A. 17b) was 
diagonalized by means of a Holstein-Primakoff transformation in Refs. [A2J. We employ here 
a closely related approach developed in [A3], which is more convenient for a generalization 
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beyond the two- level approximation and was also used in the derivation of H g D- We drop C, 
set the energy of the atomic ground states to zero, introduce the operators 



,t _ 



ek)(9k 



where qj = 2n(j /N) and j G {0, 1, . . . , N — 1}, and obtain for N — >• oo 



iV 



?^ = uxJa + ^ a^a/e + A(frJ o + 6 9o ) (a 1 " + a) 



fc=i 

N-l 



(A.18) 

(A.19a) 
(A.19b) 



In the limit of dilute excitations, that is applicable as long as the excitation energies of the 
system are finite, the b q . obey bosonic commutation relations. Note that only the j = 
collective mode couples to the radiation field. The j ^ modes are 'dark' and will be omitted 
in the following. We write b instead of b qo and arrive at 



= tia^a + Sltfb + A(6 f + b)(a) + a), 



(A.20) 



which corresponds to H g B (Eqs. (11) and ( A.41| )) with ji = 1. Later we will derive a char- 
acteristic equation for the eigenfrequencies of H g D (Eqs. ( [l2| and ( A.45b[ )). For /i = 1 this 
equation has the solutions 



24 = cj 2 + 4kxj + ft 2 ± VV 2 + 4«u; - ft 2 ) 2 + 16A 2 c^. 



(A.21) 



Diagonalization of Htm- Now we consider Htm (Eq. Q). The coupling of the electromag- 
netic field and a single harmonic oscillator 'atom' is described by 

o _v Mpi-eA) 2 mQ 2 xf 

1=1 

Note that we drop the index k numbering the atoms in % tm for a moment. As usual, we 
assume A(r) « A = Ao(a^ + a) in the region where the atoms are located. It is convenient to 
make the canonical transformation x\ = —pi/(mQ) and ^ = mQxi. This yields 

= E + ^) + ^oO(at + a) £ *, + + a)\ (A.23) 

where we have written and instead of 2^ and to keep notation simple. Succes- 
sively introducing relative and center-of-mass coordinates, {#i,.Pi, #2,^2} -Xi? Pi}, 
{Xi,Pi,x 3 ,p 3 } {£2,^2,^2,^2}, • • leads to 

<* = £(|^ + ^) + 2M + + + a)nX + ^-0 ( fl t + a) 2 . 

i=l x 7 

(A.24) 

Here, X = X n = ^ Y^=i x j an d P — Pn — Ej=i Pj are tne center-of-mass coordinates of all 
particles in the harmonic oscillator atom, and M = nm. The relative coordinates are given 

by Xi = 0-/iYl)=i x j) ~ x i+i and ^ = l/(i + l)(Yl)=iPj -iPj+i)i and V>% = ?W(i + l). Note 
that the electromagnetic field couples only to the center of mass. With this preliminary work 
done, one can write the full Hamiltonian as 

N N iV(n-l) 

H tm =cDa t a + ^E4c fe +7E(4 + c fc )(a t +a) + ^ ^ + -) + C". (A.25) 

fc=l A;=l 2=1 
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The operator c\ excites the center-of-mass degree of freedom of the fcth atom, and the gener- 
ators for the N(n — 1) relative coordinates are denoted by b\. We have introduced 




2 A 2 



>y = eAo\—, k = nN— — -, (A.26) 



and removed the ft-term by means of Co = \Juj 2 + Akuj and 7 = y/uj/ujj as before (C = (a; — 
uj + NQ)/2). The first three terms are again a special case of H g D with Qi = Q and A/ = 7 for 
all /, and \i — N. Hence, their eigenvalues follow from the roots of the characteristic equation 
for H g D (Eq. ( A.45b| )), simplified by the present conditions. They can be explicitly calculated 



and are the frequencies of the normal modes of field and center-of-mass coordinates. We find 
TV — 1 eigenfrequencies being equal to Q and represent the generators of the corresponding 
collective excitations also by b\. Only two eigenfrequencies e± are nondegenerate, 

24 = ^ 2 + 4kxj + ft 2 ± VV 2 + 4kxj - Q 2 ) 2 + 167V7 2 cj^ (A.27a) 
= cj 2 + + ft 2 ± v 7 ^ 2 + 4kxj - ft 2 ) 2 + 16A 2 cj^. (A.27b) 



We have defined A = v A/7. Since the dipole moment d of the transition from the ground state 
of an atom to its first excited state is given by d = (n\ex\n — 1) = e^n/2mQ, we can rewrite 
A = AoQd^N and n = A 2 /Q. Denoting the generators of the e± -modes by a±, we arrive at 

nN — l 

H tm = £± (aia± + -) + J2 + 2 } " I' (A ' 28) 

1=1 



Shift of the resonator frequency due to the pA- and A 2 -terms. Consider a 
system of N mutually noninteracting objects (e.g. atoms) with Hamiltonians 

^^E^f+Urf.-.i) (A.29) 

coupled to a field mode of frequency u. It is described by 

N 

H mic = coa^a + ^ (ff* ic + U k pA +n%), (A.30) 

k=l 

where ff* ic = £^ (Pi)V^i + Hnt(r{, . . . , r*J, «* A = - E?=i 9? Apf/mf, and = 
Slii fe fc ) 2 A 2 /2ra|. We denote the eigenspectrum of i^ ic by {^ fc ,|mfc}fe} and the pho- 
ton states by \l) and calculate the shifts Su p a and Sou a 2 °f the resonator frequency due to 
an d XI ^a 2 usm g the fi rs ^ nonzero terms in a perturbation series for the energy 
of |0,...,0,/). We take u <C (E^ k - E%) =: fi^ fc for m k ^ and A(r^) « A. With 
d^ fc)0 = fc(rafc| X^=i ^ r / c |0)fc? we find f° r the jth terms AE J pA and Ai£^ 2 in the perturbation 
series for the perturbations ^ and ^ 

A££ A = (A.31a) 



iV 



A^ A = _A^ ^ ^Je-d^ | + | (A .31b) 



N 



^EE^J-Col 2 (^l)-(jf)+(2I + l)(^) (A.31c) 



AE\ 2 =Al(2l + l)J2J2^k ( A - 31d ) 



fc =i i=i 2m * 
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FIG. A.l: Situation of the generalized no go theorem. Atomic spectra are drawn black, eigenenergies 
of the free electromagnetic field blue, (a) Structure of a low-energy state |£o) of the uncoupled system. 
For N — >• oo, the numbers of excited atoms £ and of photons % in \So) are small compared to N, £ <C N 
and x ^ N. (b) Structure of a component of Si. The coupling has induced one atomic transition and 
created or annihilated one photon (shown is an excitation of the second atom and the creation of a 
photon). The state \Si) is the sum of all such states. Their amplitude in the eigenstate of the coupled 
system is smaller than the amplitude of \So) by a factor oc Ao oc V -1 ^ 2 . In general, \S S ) represents 
the sum over all states obtained from \Sq) via s atomic transitions and s creations or annihilations of 
a photon. They contribute to the eigenstate of the coupled system by an amplitude oc Aq. 

Therefore, 

,2 



N 

5u pA = -2AlY J U nt k \e-dt k ^(l+-^) (A.32a) 

fe=lm fc /0 ^ L m k ) 
N n k ( n k\2 

The pA-terms cause a negative and the A 2 -terms a positive frequency shift. Note that 5uu p a 
and Suoa 2 almost cancel due to the TRK (applied for each k). The resulting total frequency 
shift Suj = SujpA + Sujaz is suppresed by ~ (gj/QJ^ ) 2 as compared with 5uj p a and 5uja 2 - 

The generalized Dicke Hamiltonian TigD- In this section, we derive the Hamiltonian 



H g D (Eq. (11)) from H m i c (in the form of Eq. (|TQ|)) for N — >• oo and show how to obtain and 
evaluate its characteristic equation. 

According to our strategy formulated in the main text, we start from low atomic densities 
and expand the eigenstates \S) of H m i c in powers of Aq oc V~ x I 2 , 

CO 

\£)ctJ2 A o\Q, (A.33) 

6 = 

where \S S ) stands for a sum over components that each describe s transitions from \Sq) both in 



its atomic and in its photonic part and hence has weight ~ Aq (Fig. A.l). The corresponding 
eigenenergies can be written as S oc ^2 SS , Ap~ s (S s \H m i C \£ s '}. We are interested only in the 
low-energy spectrum of H m ic Thus, we assume that the number of atomic excitations £ = 
J2k ^i>o \(h\£o)\ 2 and the number of photons \ — (£o|^a|£o) i n the uncoupled eigenstates 
\£q) are small compared to N, £ <C N and x <^ N. We now calculate S by dropping all 
s + s' > 2 terms and show that for the low-energy spectrum of H m i c all matrix elements that 
induce transitions in-between excited atomic states are irrelevant. To that end, we write 

N fJL nr 

Hmic = ^a + J2 (^^.^UfcX'fcl + i>io\/^(ni' — no*,^(a + + a)l'fc>('fcl) + ( A - 34 ) 

k=l l,l'=0 

with cu = \Juj 2 + 4kxj and C — (Co — cj)/2, and we define 

N fJL 

H = (ja* a + £ n i\ l k)( l k I (A.35a) 

k=l 1=1 

H cpl = (a f +a)^^AQj-Qi(ido t i\O k )(l k \ - id lt0 \l k )(0 k \) (A.35b) 

k=l 1 = 1 U 
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AH=(a^+a)^2 E A J ^ - ^(id, ^)(l k \ - id ltl ,\l k )(l' k \) . (A.35c) 



Accordingly, 

£ oc <£ \H mic \£) (A.36a) 
oc <fo|H|£b> + 4>«£b|#cpi|£i> + <fi|Hcpi|£b> + <fo|AH|fi) + <fi|Aif|£b». (A.36b) 

Let us now compare the contributions of ff cp i and AH to £. The photonic parts of H cp \ and 
AH are equal and need not be further considered. We write (H cp \) = (£ \H cp i\£i) + (£i\H cp \\£ ) 
and (AH) = (£ |Aif|£i) + (£i|AiJ|£ ), and we find 

(ah) = nwni - ^om^ gjLi ((ww) + ww))] 

(^cpi) Er=i^ Im K,iEf=i ((fb|Ofc)</fc|fi) + <5i|o fc >(z fc |£b»] ' ' 

Since A 7 ' — >• oo, the number of nonzero terms in the /c-sums is decisive. For given the 
sum over k in the numerator has at most £ nonzero terms, whereas the sum over /c in the 
denominator has at least N — £ nonzero terms. Hence, we drop Ai7, which represents the 
matrix elements of H m i c connecting the excited states of an atom, and keep H cp \ as the 
relevant coupling part of H m ic- We reintroduce the ft-term and call the resulting Hamiltonian 
generalized Dicke Hamiltonian H g r>, 

N fJL 

H gB = ooa^a + k(o) + a) 2 + ^ ^ (^|^)(^| - ^o(a f + a) (^,o|^)(O fc | + H.c.)) . (A.38) 

k=l 1=1 

It has the same low-energy spectrum as H m i c . Paralleling our treatment of Hd, we introduce 

1 N 

4,1 = \h)(0kl h\ a = —Y^ h \h){^ (A.39) 

whe re qj = 2ir(j/N) and j G {0, 1, . . . , TV - 1} as before. With ^f =1 a^a fc> / = ^ A^ z ' 

Eq. (A.38) becomes 

H gD = a;a f a + k(<J + a) 2 + ^ (fi, ^ ^.,/^-,/ - A^ l ^N{o) + a)(id Zi0 6j 0>i + H.c.)). 

z=i j=o 

(A.40) 



The operators are bosonic in the limit of dilute excitations (£ <C iV)|A3j. Again, we 
drop the energy of the dark modes, write bi instead of b q0i i, define A = AoQi\doi\y/N , and 
remove the ft-term by substituting uo — >• a; = Vcj 2 + 4kxj and A/ — >• Aj = ^uj/Co\i and adding 
C=(u>- u)/2. This gives 

/i /i 

H gD = cDa+a + ^ ftj&Jft, + ^ X^b] + ^)(a f + a) + C. (A.41) 
In order to find the eigenfrequencies of H g D, we introduce canonical coordinates by means of 
*=^(a f +a), p x =i^(a! -a), yi = -^(b\ + bi), Pi=i^Y(b]-bi), (A.42) 
and define X T = (#,2/1, . . . P T = (p x ,Pi, • • • 5 P/z), and g t = 2A|V^A- This yields 

n gD = ^ + ix T n 2 x - 1 ( w + ( A - 43 ) 
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where 



n 2 = 



(u 2 91 ■■■ g»\ 
91 n? 



(A.44) 



The orthogonal matrix G that diagonalizes induces a point transformation to the normal 
modes X = GX and P = GP. The eigenvalues ef of Q 2 are the squared eigenfrequencies of 
the system. They solve the characteristc equation 



o=(n^- e2 ))((- 2 - e2 )-E 



9t 



i=i -I 



(nw-.'))^-,')-^^) 

=1 « 



Z'=l 



None of them can be zero since this would imply 



(A.45a) 
(A.45b) 



2=1 



2m j 



(A.46) 



We have used Co = \/<^ 2 + 4/oJ, Xi = ^Juj/CjX^ \i = Ao^z |<^oz | \/^V, and k = NAq Y17=i Qi l^ m i- 
However, the left side of Eq. (A.46) is positive, whereas its right side is negative according to 
the TRK for i/° ic (Eqs. ( |A.16| ) or Eq. ^ of the main text). 



[Al] W. Thomas, Naturwissenschaften 13, 627 (1925); F. Reiche and W. Thomas, Z. Phys. 34, 510 

(1925); W. Kuhn, ibid. 33, 408 (1925). 
[A2] C. Emary and T. Brandes, Phys. Rev. Lett. 90, 044101 (2003); C. Emary and T. Brandes, Phys. 

Rev. E 67, 066203 (2003). 
[A3] J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 



